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Abstract. We present optical V and i-band light curves of the gravitationally lensed BAL quasar HE 2149—2745. 
The data, obtained with the 1.5m Danish Telescope (ESO-La Silla) between October 1998 and December 2000, are 
the first from a long-term project aimed at monitoring selected lensed quasars in the Southern Hemisphere. A time 
delay of 103±12 days is determined from the light curves. In addition, VLT/FORS1 spectra of HE 2149—2745 are 
deconvolved in order to obtain the spectrum of the faint lensing galaxy, free of any contamination by the bright 
nearby two quasar images. By cross-correlating the spectrum with galaxy-templates we obtain a tentative redshift 
estimate of z = 0.495 ± 0.01. Adopting this redshift, a P = 0.3, A = 0.7 cosmology, and a chosen analytical lens 
model, our time-delay measurement yields a Hubble constant of Ho = 66 ± 8 km s~ 4 Mpc^ 1 (la error) with an 
estimated systematic error of ±3 km s _1 Mpc -1 . Using non-parametric models yields Ho = 65 ±8 km s -1 Mpc~ x 
(la error) and confirms that the lens exhibits a very dense/concentrated mass profile. Finally, we note, as in other 
cases, that the flux ratio between the two quasar components is wavelength dependent. While the flux ratio in the 
broad emission lines - equal to 3.7 - remains constant with wavelength, the continuum of the brighter component 
is bluer. Although the data do not rule out extinction of one quasar image relative to the other as a possible 
explanation, the effect could also be produced by differential microlensing by stars in the lensing galaxy. 
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1. Introduction 


The time-delay between the gravitationally lensed images 
of a distant source is a measurable parameter. Observed 
as the time difference between the arrival dates of a single 
(lensed) wavefront emitted by a distant source, it is di¬ 
rectly related to the Hubble constant H 0 (Refsdal 1964). 
Obtaining accurate time-delay measurements in multi¬ 
ply lensed quasars can therefore yield (i) a determina¬ 
tion of Ho provided the mass distribution in the lens 
is known, or (ii) constraints on the mass distribution 
in a given lens, using Hq as inferred from other meth¬ 
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ods. During the last 20 years much effort has been de¬ 
voted to the observations of lensed quasars, and in par¬ 
ticular to long-term monitoring of selected systems. Some 


of these are Q0957+561 (Schild 199C, Vanderriest et al. 
|l989|) , PG1115+0 80 (Sc hechter et al. |l997|) , B1608+656 
(Fassnacht et al. 1999| ) and B1600+434 (Burud et al. 
|2000| and Koopmans et al. 2000| ). In this context, we 
have been conducting a photometric monitoring pro¬ 
gram at the Danish 1.5-m telescope at La Silla ob¬ 
servatory (ESO, Chile) since October 1998, with the 
goal of measuring the time-delays in several well studied 
lensed quasars. We present here the first result from this 
program: the time-delay measurement in the two-image 
quasar HE 2149—2745. 

The lensed nature of the BAL quasar HE 2149—2745 
at z = 2.03 was established by Wisotzki et al. ( 1996 ). 
This system proves to be an easy target for monitoring 
at a site with reasonable seeing conditions (up to 2"): it 
is bright (B = 17.3) and the two quasar images have an 
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Fig. 1. Field of view 5'x5' in size around HE 2149—2745. 
The three reference stars (labelled SI, S2, S3) used for the 
photometry and the two PSF-stars (PSF1, PSF2) used in 
the spectral deconvolution are indicated. North is up and 
East to the left. 

Fig. 2. Left: Stacked F-band image of a 9"x 9"region cen¬ 
tered on HE 2149—2745. The seeing is ~2" and the total 
exposure time is ~ 17 hours. Right: Deconvolved image 
(FWHM=074) obtained from the simultaneous deconvo¬ 
lution of 57 frames. North is up and East is to the left. 


angular separation of 1.7". The monitoring program, the 
light curves and the time delay are discussed in Sect. I- 0 
and |] below. 

The lensing galaxy has been detected in HST NICMOS 
and WFPC2 images but its redshift remains unknown. 
With the aim of measuring this redshift, we have obtained 
a spectrum of HE 2149-2745 with FORS1 at UT1 (ESO- 
Paranal, Chile). The analysis of the spectroscopic data is 
described in Sect. |[ This section also includes a discussion 
on the spectral differences between the two quasar compo¬ 
nents. Mass models and estimates of the Hubble constant 
are presented in Sect. [| Finally, Sect. [?] summarises the 
main results. 


2. Photometric Monitoring at the 1.5m Danish 
Telescope 

2.1. Observations and Data Reduction 

Weekly observations of HE 2149—2745 were carried out at 
the Danish 1.5-m telescope at ESO-La Silla from October 
1998 to December 2000. The target is visible from the 
beginning of June to the end of December, which results 
in gaps of about ~ 5 months in the light curves. Apart 
from these gaps, very few points were missed because of 
poor weather or technical problems. Observations were ob¬ 
tained in the F and Gunn i bands with DFOSC (Danish 
Faint Object Spectrograph Camera) which has a pixel size 
of 07395. Fig. [l] shows the central region of the held of view 
of DFOSC, i.e., 5 arcmin 2 of the 13 arcmin 2 available. The 
F-band was given priority because of the detector’s bet¬ 
ter sensitivity at these wavelengths, but observations were 
also carried out in the j-band in order to monitor possible 
colour changes. The exposure time was set to 900 sec dur¬ 
ing the first 6 months of observations. It was increased to 
1800 sec for the rest of the observations in order to improve 
the photometric quality for the faint B quasar image. For 
each observed data point the total exposure time was di¬ 
vided into three dithered exposures in order to remove 
cosmic rays and bad pixels. The seeing varied from 077 to 
279, with a median of 176. All the imaging data were pre- 
processed (bias-subtracted and flat-field corrected) using 
standard IRAF routines. 


Table 1. Positions of three reference stars relative to the 
position of the A component in HE 2149—2745. 


Star 

R.A. (arcsec) 

DEC (arcsec) 

SI 

-48.126 

136.268 

S2 

-57.126 

-72.647 

S3 

-38.160 

-4.145 


3. Photometry 


The photometry of the blended quasar images was per¬ 
formed by applying the MCS deconvolution algorithm 
(Magain, Courbin & Sohy 1998 ). This algorithm has al¬ 
ready been used to analy ze the data of seve ral le nsed 
quasars (e.g., Burud et al. 2000 , Hjorth et al. 2001 ). Its 
main advantage is its ability to use all the data, even 
rather poor, irrespective of image quality and lunar phase. 
The final deconvolved image is produced by simultane¬ 
ously deconvolving all the individual frames of the same 
object from all epochs. The positions of the quasar im¬ 
ages (two in the present case) and the shape of the lensing 
galaxy are the same for all the images and are therefore 
constrained using the total S/N of the whole data-set. The 
intensity of the point sources are allowed to vary from im¬ 
age to image, hence producing the light curves. 

The two quasar components are well separated in our 
deconvolved image of HE 2149—2745 (see Fig. ||), but the 
lensing galaxy remains too faint to be detected. This is for¬ 
tuitous since contamination by an extended object would 
complicate the analysis. 

The light curves of HE 2149—2745 consist of 57 data 
points in the F-band and 41 points in the Gunn i band, as 
presented in Fig. The points between JD 2451800 and 
.ID 2451900 were obtained just after the installation of the 
new chip on DFOSC. The new chip does not have prob¬ 
lems of charge diffusion, as the old chip had, and conse¬ 
quently, the effective seeing is now considerably improved. 
Our original PSF star, for which the ADU counts were 
always in the linear regime on the old CCD, became sat¬ 
urated on the new chip. Another PSF star further away 
from the target had to be used on a few frames before the 
exposure time was adjusted to the new chip. These points 
are therefore more noisy than the previous points. 

The data are plotted relative to a reference star com¬ 
mon to all images, star S2 in Table [j]. Two other stars were 
deconvolved as well, in order to check the relative photom¬ 
etry and to check for systematic errors (see Table [I] for the 
positions of the stars). The error-bars include both photon 
noise and additional systematic errors, e.g., PSF errors. 
The latter is estimated by using one of the reference stars, 
as explained in Burud et al. (2000). 


4. Time-delay measurement 

By sliding the light curves across one another, one can 
make a rough ”by eye” estimate of the time-delay of ~ 100 
days, with A leading B. Usi ng the % 2 minimization method 
described in Burud et al. ( 2000 ), a more objective value 
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Fig. 3. Light curves of the two images A and B of HE 2149—2745 and the comparison stars SI and S3 (Table |lj). 
The F-band light curves are displayed on the left while the Gunn z-band light curves are shown on the right. The 
magnitudes of the quasar and stars are calculated relative to the reference star S2. The error bars include photon noise 
and PSF errors estimated from the deconvolution of a reference star. For display purposes, the B component is shifted 
by —1.3 mag. and —1.2 mag. in the F and Gunn i band, respectively, with respect to their original values. Likewise, 
SI is shifted by 1.1 and 1.35, and S3 is shifted by +1.15 and +1.35 mag in the F and Gunn i filters, respectively. S3 
was saturated in the last Gunn i images and is therefore not shown on the plot for these frames. 


of At = 103 ± 12 days is found from the F-band light 
curves. The z-band data are noisier and have fewer data 
points but the measured time-delay, At = 104 ± 31 days, 
agrees well with the delay derived from the F-band. A 
simultaneous minimization of the F and i band curves 
gave a A t = 109 ± 22 days. The errors quoted here are 
obtained from Monte Carlo simulations of 1000 sets of 
light curves, assuming that the photometric errors are 
uncorrelated and follow a Gaussian distribution. The y 2 
minimization method was also performed on 1000 sets of 
curves where 5 randomly chosen data points were removed 
from each set. The results from these simulations yielded 
a time-delay of At = 101 ± 30 days and At = 107 ± 30 
days in F and i respectively, confirming that the time- 
delay measurement is robust. The magnitude differences 
between the A and the B components are found to be 
1.635 ± 0.001 mag in F and 1.505 ± 0.003 mag in the i- 
band. This corresponds to flux ratios of 4.51 and 4.00 in F 
and i bands respectively Although no erratic changes are 
detected in the light curves within the measurement er¬ 
rors, microlensing on long time scales may still be present 
as suggested by the VLT spectra presented in Sect. ||). 
Given that the F-band data contain more points and are 
less noisy than the z-band data, we adopt the F-band es¬ 
timate of the time-delay as the best one: At = 103 ± 12 
days (1(7 error). Fig. || summarizes the photometric data 
and Fig. ^ displays the F and z-band light curves with 
component B shifted by 103 days. 


5. Spectroscopy of HE 2149-2745 

5.1. VLT spectroscopy 

To convert the time-delay into an estimate of the Hubble 
constant, one must know the geometry of the lensing 
system and this includes the redshift of the lensing 
galaxy. For this purpose, we took an optical spectrum 
of HE 2149-2745 with the ESO VLT/UT1 using the 
Multi-Object-Spectroscopy (MOS) capability of FORS1. 
The observations were obtained on November 19, 2000 
under fairly good seeing conditions (0.8"). The 1" slitlets 
of FORS1 were aligned to obtain simultaneously the spec¬ 
trum of HE 2149—2745 and of 2 PSF stars (indicated in 
Fig. |) about as bright as the two quasar images. The high 
resolution collimator was used, resulting in a pixel size of 
0.1", in combination with the G300V grism and GG435 
order sorting filter. Three exposures, each of 1000s, were 
taken so that cosmic rays could be removed. 

5.2. Redshift of the lensing galaxy 

Using the spatial information in the spectra of the two 
PSF stars, the spectrum was spatially deconvolved with 
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JD—2451000 


Fig. 4. The V (top) and Tband (bottom) light curves where the B curve is shifted by 103 days. The light curve of 
component B is shifted by 1.635 and 1.505 magnitudes in the V and i-band, respectively, with respect to their original 
values. The curves for the A component are as in Fig. 


the spectral version of the MCS deconvolution method 
(Courbin et al. {2000a ). The deconvolution process de¬ 
composes the data into the individual spectra of the two 
quasar images and the faint lensing galaxy (see Figs. || 
and |7|), in a way similar to the one described in Lidman 
et al. fcoQOD for HE 1104-1805. 


The lens spectrum was cross-correlated with a tem¬ 
plate spectrum of an elliptical galaxy (Kinney et al. 1996| ). 
Due to the low S/N the calculated correlations are low, 
nevertheless, there is a significant “bump” in the correla¬ 
tion function covering the redshift range 0.49 < z < 0.60 
as shown in Fig. The spectrum of the lens and the tem¬ 
plate are shown in Fig. ||[ where the template has been 
shifted vertically for clarity. The Calcium H and K lines, 
the G-band and the Mg triplet are labeled. A measure 
of the reliability of the redshif t esti mate is given by the 
r-statistic of Tonry and Davis ( 1979 ), which is the signal- 
to-noise ratio of the main peak in the cross correlation. 
We find that r = 1.9 meaning that the signal-to-noise of 
the correlation peak is poor. A value above 3 is considered 
secure (Kurtz & Mink 


As can be seen on the HST image (cf. Fig. 0 ) the 
lensing galaxy lies very close to the B image. We therefore 
looked for evidence of absorption in the spectrum of the 
background quasar. In order to remove the quasar spec¬ 
tral features, the B spectrum was divided by that of the A 
component (Fig. ||). Two, possible Ca H & K candidates 
were found, one at z = 0.489 and another at z = 0.504. An 
integrated A+B spectrum at similar resolution but cover¬ 
ing sh orter wavelengths has been published by Wisotzki 
et al. 1996 We looked for possible Mgll absorption in 


the quasar spectrum at the putative redshift of the lens. 
Unfortunately, the 2 = 0.489 line would fall in a strong ab¬ 
sorption line of the quasar, while a candidate at z = 0.504 
is found in the wing of a quasar emission. The signal-to- 
noise ratio of the lens galaxy and the Ca H and K ab¬ 
sorption features in the B component are very low and we 
can not rule out any of the two estimates. We therefore 
take z = 0.495 ± 0.01 to be the most likely redshift for the 
lensing galaxy. This is within the estimated redshift range 
0.37 < z < 0.50 which is based on the position of the lens 


in the fundamental plane (Kochanek et al. 2000). 


1998|) . 















































Burud et al.: An Optical Time-Delay for the Lensed BAL Quasar HE 2149—2745 


5 



Fig. 5. Extracted spectrum of the lensing gala xy (bo 
and the template spectrum from Kinney et al. ( 199G| ) 
used in the cross correlation. The main spectral fea 
Calcium H and K lines, the G-band and the Mg ti 
are indicated by dashed vertical lines. The lens reds 
tentatively z = 0.489. 



Fig. 6. Cross-correlation function of the galax y spe ctrum 
with the galaxy template from Kinney et al. ( 1996 ). The 
centre of the correlation peak is at z = 0.535. 



Fig. 7. Top: spectra of the A and B components. Bottom: 
the spectrum showing A—3.7xB. The scaling factor is cho¬ 
sen to best cancel the emission lines. Telluric absorption 
lines are indicated. 



Wavelength (A) 


Fig. 8. Spectrum of the B component divided by that of 
the A component of HE 2149—2745. The plot shows a 
zoom on the region corresponding to possible Calcium H 
and K absorption due to the lens at z = 0.489. The bump 
between 5750 and 5900 A is the residual CIIIJ+A1III emis¬ 
sion from the quasar. 


5.3. Spectroscopy of the quasar images: colour 
differences, extinction or microlensing ? 

The quasar spectra show that the continuum of com¬ 
ponent A is bluer than the continuum of component B 
(Fig. |). This confirms the difference in flux ratio that 


was found from the light curves in the V and Gunn i- 
bands (Sect. [|). We note three points of interest about 
these spectral differences: 
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JD—2451000 

Fig. 9. Colour excess of the B component relative to the 
A component, (V — I)b — {V — I) a as a function of time. 
No significant variations can be detected over a period of 
900 days. The solid line shows the best fit line through the 
points. 

1. As shown in Fig. the colour difference between the 
quasar images does not vary during the period of our 
observations (~ 3 years). 

2. The emission line flux ratios at the date of the ob¬ 
servation are not wavelength dependent: all emission 
lines cancel in the difference spectrum of Fig. This 
difference spectrum has been produced by multiplying 
the spectrum of the B component by a constant fac¬ 
tor of 3.7 and by subtracting it from the spectrum of 
the A component. Fig. [To| gives a different representa¬ 
tion of the same effect, where the spectra of the two 
quasar images are normalized to the same continuum. 
The equivalent widths of all emission lines (including 
the Fell pseudo-continuum) are smaller in A. 

3. The CIV broad absorption line behaves as the contin¬ 
uum: its equivalent width is unchanged in the A and 
B spectra (Fig. |Tcj ). 

The spectral differences cited above have essentially 
two possible explanations: (i) differential magnification 
due to microlensing, or (it) differential reddening of the 
B component by the lensing galaxy. 


5.3.1. Microlensing 

If microlensing is present in HE 2149—2745, one expects 
the continuum region and the much more extended emis¬ 
sion lines re gion t o be affected differently (Wambsganss 
& Paczynski 1991 ). In most AGN models, the size of the 
continuum region depends on wavelength, and the Broad 
Emission Line Region (BELR) is usually believed to be 



Wavelength (A) 

Fig. 10. The A (thick) and B (thin) spectra of 
HE 2149—2745 normalized to a common continuum. The 
equivalent widths of the emission lines, including the Fell 
pseudo-continuum, are smaller in A. The equivalent width 
of the CIV broad absorption line (~4650 A) is the same 
in both the A and B components. 


larger than the continuum region by more than one or¬ 
der of magnitude. Since the B component is closest to the 
lensing galaxy we would expect that microlensing effects 
occur most often in the B-conrponent. Microlensing of the 
A component however is also possible. 

Assuming the lensed quasar in HE 2149—2745 follows 
the “standard” AGN model and that we are observing 
component A through a network of caustics produced by 
stars in the main lensing galaxy, one can imagine a sce¬ 
nario where the inner - and bluer - continuum region of 
the component A is being enhanced by a larger amount 
than the outer redder parts. As the BELR is much larger 
than the central AGN (even the redder parts), it remains 
unaffected by microlensing. This interpretation has al¬ 
ready been proposed to explain similar spectral differences 
observed in the doub le HE1 104-1805 (Wisotzki et al. 1993, 
1995, Courbin et al. 2000b ). 

In the case of HE 2149—2745 the Einstein radius of 
a typical deflector in the lens plane projected onto the 
source plane is (H = 0.3, A = 0): 


Re = 


AGM D ds D a 


D d 


= 1.3 x 1O- 3 (M/O.1 M 0 ) 1/2 P c (1) 


where D d , D s and D ds are angular diameter distances to 
the deflector, the source and between deflector and source 
respectively, and M is the mean mass of all microlenses. 
The duration of crossing the Einstein radius will be up 
to 10 years assuming a relative velocity between source 
and lens that is equal or larger than the velocity in the 
lensing galax y whi ch is of the order of ~ (100 — 400)kms _1 
(Mould et al. 1993| ). Although shorter microlensing events, 
such as a caustic crossings, can occur, the long time scale 
of crossing the Einstein radius shows that an event may 
be present in our data, and stable over our relatively short 
period of observation of 900 days, as Fig. p] would suggest. 
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Fig. 11. The residual magnification [i of A relative to B 
as a function of wavelength. This magnification factor is 
derived by fitting the A/B continuum ratio after correct¬ 
ing the B spectrum by the constant factor 3.7 (as in the 
lower panel of Fig. [5]) 



Fig. 12. The magnification of component A in magnitude 
(Am = 2.5 log//) as a function of inverse wavelength (solid 
line). The wavelength dependence is well reproduced by 
Am = 0.145 ) (dashed line) 


Assuming that the continuum of A is amplified by mi- 
crolensing, we may derive the wavelength dependence of 
the magnification from the A/B continuum ratio. This 
relative magnification factor /i is given in Fig. |ll| It is de¬ 
rived by fitting the A/B continuum ratio after correcting 
the spectrum of B by the constant factor 3.7. The relative 
magnification is directly related to the size of the contin¬ 
uum region as a function of wavelength (e.g. Schneider 
et al. 1992), the precise measurement of which requires 
the knowledge of the background magnification fiQ due to 
other distant subimages. While /i can be estimated from 
Fig. 0 fi ,o can only be obtained by measuring the back¬ 
ground magnification before and after a caustic crossing 
event. This can be done with a spectrophotometric moni¬ 
toring. In Fig. 0 we have plotted the magnification fac¬ 
tor - expressed as a magnitude difference - as a function 
of the inverse wavelength. It shows a remarkably tight 
wavelength dependence strikingly similar to that found by 
Nadeau et al. (1999) for Q2237+0305. This suggests that 
HE 2149—2745 presently suffers a chromatic microlensing 
event similar (although of much longer duration) to the 
1991 high-magnification event in Q2237+0305 and which 
was interpreted - on the basis of its colour dependence - 
as evidence for the thermal accretion disk model as the 
source of UV-visible continuum emission in quasars (cf. 
Nadeau et al. 1999). 



Wavelength (A) 


Fig. 13. An enlargement of the spectra of components A 
(upper) and B (middle, multiplied by 3.7), together with 
the difference spectrum A—3.7xB (lower, multiplied by 
2). The region of the CIII] emission line is illustrated. The 
wavelengths of A1III AA1855, 1863 and CIII] A1909 are 
indicated at the redshift of the quasar (z= 2.03). This 
blend is typical of BAL quasars (Hartig & Baldwin 1986). 
Some subtle differences in the blue part of the line profiles 
may be noticed (see text). 


While broad emission lines should be essentially un¬ 
affected by the caustic crossing, some subtle differences 
may arise in the line profiles if some parts of the BELR 
are selectively magnified (Nemiroff 1988, Schneider & 
Wambsganss 1990). Looking in detail at the line profiles 
of components A and B (Fig. 0), small differences are 
indeed present in the blue wing of the A1III and CIII] 
emission lines. Interpreted in terms of microlensing, this 
could indicate that part of the BELR is being magnified. 
Similar profile variations are predicted by the models of 
Schneider & Wambsganss (1990). Interestingly, no differ¬ 
ence in the CIV BAL profile, like those reported for the 


quadruple BAL quasar H1413+117; cf. Angonin et al. 1990 
and Hutsemekers 1993, is noticed. 

The present findings may be of importance for investi¬ 
gating the inner structure of quasars. Spectroscopic data 
obtained at regular time intervals, e.g., spectra separated 
by intervals corresponding to the time-delay, will test the 
microlensing interpretation, cancel possible time-delay ef¬ 
fects (see Wisotzki et al. 1995) and help to probe the 
innermost region of quasar structure, using microlensing 
magnification. 


5.3.2. Reddening 

An alternative explanation for the observed colour differ¬ 
ence between quasar images is that the B image is red¬ 
dened by the lensing galaxy. Assuming this is true, we 
can estimate the amount of reddening which would corre¬ 
spond to the strength of the putative Ca H & K absorption 
lines seen in the B spectrum (Fig. |j), if average Galactic 
conditions apply to the present case. From the interstel¬ 
lar line observations of Sembach et al. (1993), we derive 
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a mean relation between E(B-V) and the Ca K equiva¬ 
lent width. With a de-redshifted equivalent width ~ 100 
mA, this would imply a colour excess E(B-V) ~ 0.05, 
which corresponds to E(V-I) ~ 0.1, using the extinction 
curves of Mathis (1990), after proper redshift corrections 
are applied. This is compatible with the measured colour 
difference of 0.15 (Fig. |j). The extinction can also be esti¬ 
mated by fitting the magnitude difference of the A and B 
components as a function of wavelength (Fig. |l2|) with a 
typical extinction law. In the wavelength range of interest, 
the Galac tic and SMC extinction curves are very similar 
(Pei |l992|) and can be reasonably well represented by 


Ax = 


1.75 E{B - V) 


( 2 ) 


where A\ is the extinction in magnitudes, and where the 
ratio of the total-to-selective extinction Ry is taken to be 
Rv = 3.1. Identifying the wavelength dependence of the 
magnitude difference in Fig. 0 to this extinction curve at 
the redshift of the lens, we have 1.75 E(B-V)(l+ 2 i; e „ s ) = 
0.145 i.e., E(B-V)= 0.05 which is in agreement with what 
is found from the Ca H & K lines. 

In the case of extinction by dust in the lens we would 
expect the continuum and the BELR to be equally af¬ 
fected. This does not seem to be the case; however, the 
wavelength range between CIV and the CIIIJ-A1III com¬ 
plex is small. A more sensitive test of the dust hypothesis 
will be to observe H-alpha in the infra-red, as this will pro¬ 
vide a larger wavelength baseline to measure differences. 

On the basis of the available data, no interpretation 
of the apparent reddening can be excluded, although mi- 
crolensing seems a somewhat more natural explanation for 
the observed spectral differences. 


6. Lens modeling 


The time-delay measured for HE 2149—2745 can be used 
to infer an estimate of Ho, based on modeling of the total 
gravitational potential responsible for the lensing effect. 
This includes the main lensing galaxy and any interven¬ 
ing massive cluster along the line of sight to the quasar. 
HE 2149—2745 is a relatively easy case to model, as there 
is no obvious mass concentration along the line of sight 
to the source apart from the main lensing galaxy. Only 
a marginally detected galaxy concentration is seen 30" 
West and 1' North of the main lens. In addition, the lens¬ 
ing galaxy is almost aligned with the two quasar images 
(see Figs. |l| and 0) . Two different approaches can be used 
to convert time delays into Ho: models that involve an 
analytical form for the lensing galaxy, and models that 
involve more degrees of freedom, e.g., a pixelated surface 
mass density (Saha & Williams 1997 Williams & Saha, 
200C). We use both approaches. 


6.1. Analytical models 


Our ana lytic model uses the mass profile of Kneib et 
al. (1996): a truncated Pseudo-Isothermal-Elliptical-Mass- 



1 


2 1 0-1 

a (arcsec) 

Fig. 14. Field of 3" around HE 2149—2745, where the 
lensing galaxy is indicated as G, close to the faint quasar 
image B. Also indicated are the critical and caustic curves 
created by a single lensing galaxy and no external shear. 
The image has been obtained from the CASTLEs public 
survey of gravitational lenses. This I-band image comes 
from HST/WFPC2 (Kochanek et al.). 


Distribution (PIEMD). Truncated PIEMD are elliptical 
mass distributions smoothly truncated at radius r cu t . For 
radii smaller (respectively larger) than r cu t, the projected 
surface mass density profile is varying as r^ 1 (respectively 
r~ 3 ). The interest in such profile is their ability of charac¬ 
terizing any ellipticity in the mass distribution as well as 
having a total finite mass. We fit the model to the publicly 
available HST/WFPC2 data (which offers better resolu¬ 
tion and sampling than NICMOS) obtained by Kochanek 
et al., and use the emission line flux ratio of 3.7 calculated 
from our VLT spectra. During the fit, the astrometry of 
the quasar images relative to the lens is fixed, as well as 
the redshift of the source and lens. The free parameters 
include the velocity dispersion, ellipticity, which is here 
defined as [a 2 — b 2 ]/[a 2 + b 2 ] where a and b are the long 
and short axes of the lens, and cut radius of the lens. 

A common problem in lens modeling is the mass-sheet 
degeneracy. In other words, several mass profiles and in 
particular mass profiles of different compactness, will re¬ 
produce the same time-delay (Gorenstein et al. 1988). One 
has therefore to constrain the models using observational 
data and explore the widest range of physically acceptable 
models. A first constraint comes from the HST images: the 
lens is almost round. After PSF-subtraction of the quasar 
images we fit the shape of the lensing galaxy and obtain 
an ellipticity of about 0.1-0.2. Assuming that the mass 
distribution follows the light distribution we can restrict 
ourselves to this range of ellipticities. Fig. 0 shows a fam- 
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Fig. 15. Results of the model fitting when exploring the 
ellipticity/cut radius space. The lines show the best mod¬ 
els, i.e., they show the area where the models fit the data 
at lcr (inner lines) and 3cr (outer lines). Values for Ho are 
given for different models among the most likely. The hor¬ 
izontal lines give the plausible range of ellipticities, given 
the HST images. The plot displayed here uses Q = 0.3, 
A = 0.7. 


ily of models exploring the ellipticity vs. r cu t space and we 
see that little spread in the values of H 0 is observed in the 
range of ellipticities 0.1-0.2. We assume ellipticity e = 0.15 
which implies r cu t = 1.6" and Ho = 66 ± 8 km s _1 Mpc -1 
with additional systematic errors of ± 3 km s^ 1 Mpc -1 
due to the limited range of ellipticity. We also try to in¬ 
vestigate the existence of a galaxy core. For this we scan 
the possible values for the core radius. The lens model ex¬ 
cludes any value of r core larger than 0.02 arcsec. Within 
this limit, the core radius has basically no effect on the 
determination of Ho, so the choice of a given core radius 
is not critical. In our fiducial model, we use a core radius 
of ~ 0.002 arcsec which corresponds to ~10 pc. Using 
this, and keeping the other parameters free, we find a lens 
with the parameters summarized in Table We show in 
Table || the effect of a change in cosmology. Due to the 
relatively low redshift of the lensing galaxy the influence 
of a change in cosmology is small compared to the mea¬ 
surement errors. 

If we take the currently most popular cosmology, i.e., 
f2 = 0.3, A = 0.7, and A t = 103 ± 12 days, we find 
Ho = 66 ± 8 km s -1 Mpc -1 with an additional estimated 
systematic error of ± 3 km s^ 1 Mpc -1 due to parameter 
choices. 


Table 2. Model parameters for HE 2149—2745. The val¬ 
ues between brackets are fixed during the fit, the others 
are fitted. The core radius is negligible. 


Parameter 

Fit 

H, A 

[0.3, 0.7] 

Ellipticity 

0.13 ± 0.05 

P.A. (degrees) 

62 ± 5 

Vel. disp. (kms -1 ) 

198 ± 10 

r cut 

[1.25"] 

Core radius 

[0.002"= 0.01 kpc] 

time-delay 

[103 days] 


Table 3. Influence of the cosmology on the modeling of 
HE 2149—2745. Values of Ho are listed for an Einstein- 
deSitter Universe {il = 1.0, A = 0.), an open Universe 
(Q = 0.3, A = 0.) and a flat Universe with non-zero cosmo¬ 
logical constant (Q = 0.3, A = 0.7). The model parameters 
are fixed and set to the ones in Table |j. 


Cosmology 

Ho (km s 1 Mpc x ) 

H = 0.3, A = 0.7 

66 ± 8 

H = 0.3, A = 0.0 

67 ± 8 

Q = 1.0, A = 0.0 

61.5 ± 8 


6.2. Pixelated models 


The pixelated models o f Sa ha & Williams fll997| ) 
and Williams & Saha ( 2000| ) were also applied to 
HE 2149—2745. While the mass distribution resulting 
from these models is not as well related to the physical 
parameters of the lens, they allow us to explore a wide 
range of lens shapes and to estimate the robustness of the 
time-delay conversion. Fig. E shows the mass distribu¬ 
tion found by the models. During the fit, the ellipticity 
of the mass distribution was kept free. The program re¬ 
constructs almost round mass distributions, as does the 
analytical model (Fig. |l5|). We also examine the impact 
of a change in the compactness of the mass profile on the 
derived value for H 0 . Fig. E shows the result of a Monte- 
Carlo simulation where 100 models were run with different 
indices for the lens shape. The density index in Fig. E rep¬ 
resents the logarithmic projected density gradient in the 
vicinity of the images (cf. Williams & Saha, 2000 ). The 
upper panel of the figure shows that most models predict 
relatively steep/concentrated mass profiles. 

In the lower panel of Fig. 17 we display the Ho prob¬ 
ability distribution. The mean of the distribution peaks 
at Ho = 65 ± 15 km s^ 1 Mpc -1 at the 2-cr level and 
Hq = 65 ± 8 km s^ 1 Mpc -1 at the l-o level. 


7. Summary - conclusions 

We have presented the first result of a long-term photo¬ 
metric monitoring campaign undertaken at ESO between 
1998 and 2000 at the 1.54-m Danish telescope. Our V 
and i light curves allow us to measure a time-delay of 
At = 103 ± 12 days between the two quasar images of 
HE 2149-2745. 
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Fig. 16. Mass distribution found with the pixelated mod¬ 
els with cosmology Q = 0.3,A = 0.7. Almost no ellipticity 
is found. The contour levels represent k = 1/3, 2/3, 1, etc. 
The scale is 3" on a side and the astrometry is the same 
as for the analytical models. The dashed ellipses around 
the quasar images represent the magnification matrix, i.e., 
the axes are proportional to the eigenvalues and the orien¬ 
tation corresponds to the eigenvectors. The spot between 
the quasar images is the position of the lensed source in 
the source plan. 


From VLT spectroscopy, we have derived a tentative 
estimate of the lens redsliift to be z = 0.489. Applying 
both analytic and numerical lens models to the case of 
HE 2149—2745, we derive Ho = 66 ±6 km s -1 Mpc -1 with 
an additional systematic error of ± 3 km s -1 MpW 1 in the 
case of the analytic models, due the limited range of the el¬ 
lipticity. The derived mass models are relatively compact, 
although not as compact as the light profile of the galaxy. 
An extra source of systematics might be introduced by the 
uncertainty on the lens redshift estimate. This is however 
not critical for the determination of Ho as the error is dom¬ 
inated by the uncertainty in the time-delay measurement. 
As HE 2149—2745 shows smooth light curves, it is likely 
that the situation can be improved by continued monitor¬ 
ing. With an improved time-delay it would also be highly 
desirable to re-determine the lens redshift more precisely. 

Our monitoring program of HE 2149—2745 is the first 
to be carried out in two bands on such a regular basis and 
over such a long time scale. Given the error bars, we do not 
see any significant colour variation over the 900 days of ob¬ 
servation. Moreover, our spectra of the two quasar images 
show that the flux ratios in the broad emission lines be¬ 
have differently for the continuum flux ratio, and that the 
flux ratio measured in the BAL structure of the source fol¬ 
lows the behavior as the continuum region. Such behavior 



0 0.2 0 4 0.6 0.8 


h 



0 0.2 0.4 0.6 0.8 1 
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Fig. 17. Top: radial index profile vs. Ho for one hun¬ 
dred different non-parametric models, showing that all the 
galaxy models have a steep mass profile (< -0.5) Bottom: 
probability density for Ho for the same models. Units on 
the X axis of both figures is in ft-ioo • The median value of 
the distribution is 65 km s^ 1 Mpc -1 . The black dot marks 
the location of an isothermal mass profile (radial profile 
index = -1). 

can be explained both by microlensing or by differential 
extinction by the lensing galaxy, or both. So far, the data 
do not allow us to distinguish between the two possible 
explanations. In order to confirm microlensing, one would 
need for example to know the time-scale of the putative 
event, i.e., to measure the absolute magnification of the 
event and its duration. This would not only allow us to 
confirm microlensing, but also to use it to map the radial 
structure of the central AGN in the source. Although the 
Einstein radius crossing time is long for HE 2149—2745, 
of the order of 10 years, much shorter microlensing events, 
such as a caustic crossings, can occur. Conducting a long¬ 
term spectrophotometric monitoring could therefore allow 
us to probe the AGN size HE 2149—2745. 
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